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SUMMARY 


This  report  summarizes  work  porformed  during  the  concluding  phase  of  a  theoretical 
investigation  into  the  physics  responsible  for  the  initiation  and  propagation  of  Laser 
Absorption  Waves  that  have  been  observed  to  form  over  opaque  solid  targets  irradiated 
in  air  by  an  intense  laser  beam.  The  report  deals  separately  with  the  initiation  of  sub¬ 
sonic  Laser  Supported  Combustion  (LSC)  Waves  over  metal  targets  irradiated  by  a  CO„ 
laser  at  intensities  up  to  a  few  MW/cm  for  times  of  the  order  of  milliseconds  or  longer; 

and  with  the  "prompt"  initiation  of  supersonic  Laser  Supported  Detonation  (LSD)  Waves 

2 

by  short-pulse  irradiation  at  COg  laser  intensities  of  the  order  of  30  MW/cm  or  greater. 
In  the  former  case,  we  had  shown  previously  that  the  initiation  of  subsonic  waves  is  de¬ 
pendent  on  large-scale  vaporization  of  target  material,  laser-induced  inverse  brems- 
strahlung  heating  of  the  vapor,  and  subsequent  heating  of  ambient  air  by  thermal  reradi¬ 
ation  from  the  heated  vapor.  Recent  refinements  of  the  basic  theoretical  model  are 
presented  herein,  and  are  used  as  a  basis  for  correlating  existing  experimental  data  on 
LSC  wave  formation. 

Two  mechanisms  that  could  be  responsible  for  prompt  initiation  of  laser-supported  ab- 
sorption  waves  in  air  at  C02  laser  beam  intensities  above  about  30  MW/cm  have  been 
investigated.  These  are  high-field  electron-emission  from  protrusions  on  the  target 
surface,  and  thermionic  electron-emission  from  anomalous  patches  on  the  target  sur¬ 
face.  Calculations  of  absorption  wave  initiation  time  based  on  a  thermally  augmented 
high-field  electron-emission  model  are  not  inconsistent  with  experimental  measurements; 
however,  the  general  trend  of  the  experimental  data  is  not  well  predicted  at  the  lower 
end  of  the  intensity  range  near  threshold.  Calculations  of  absorption  wave  initiation  time 
based  on  the  thermionic  electron-emission  model,  with  suitably  chosen  parameters, 
agree  quite  well  with  available  exper'mental  data.  Either  of  two  widely  differing  assump¬ 
tions  about  the  detailed  behavior  of  the  intensely-irradiated  target  surface  can  lead  to  the 
aforementioned  thermionic  predictions.  One  is  the  presence  of  surface  lamina,  or 
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flakes,  that  are  thermally  isolated  from  the  underlying  material  and  that  are  of  the 
order  of  the  electromagnetic  skin  depth  in  thickness.  The  other  assumption  is  that  heat, 
conduction  into  thp  bulk  of  the  target  material  is  rate-limited,  so  that  the  temperature 
rise  due  to  laser  power  absorption  within  the  skin-depth  layer  at  very  high  incident 
power  densities  is  affected  only  by  the  thermal  capacity  of  the  material.  This  is  a 
heretofore  uninvestigated  phenomenon,  which  has  far-reaching  fundamental  implications. 
Under  this  assumption,  the  thermionic  electron-emission  initiation  mechanism  involves 
no  additional  assumptions  nor  arbitrarily  adjustable  parameters,  but  involves  only  the 
intrinsic  optical  and  physical  properties  of  the  target  material. 

In  addition  to  the  described  analyses  of  laser  Absorption  Wave  Formation  physics,  we 
have  continued  our  earlier  formulation  of  an  analysis  of  a  two-dimensional  subsonic  LSC 
wave  propagating  along  a  finite -diameter  laser  beam.  The  formulation  is  used  as  the 
basis  for  an  approximate  solution  that  shows  the  significant  effect  of  radial  flow  on  the 
wave  propagation  speed. 
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Soctlon  1 

INTRODUCTION  AND  SUMMARY  OF  PREVIOUS  WORK 


The  effort  described  in  this  report  and  in  a  previously  published  Midterm  Report 
(Ref.  1)  represents  the  continuation  of  an  earlier  theoretical  investigation  (Refs.  2 
and  3)  into  the  physics  responsible  for  the  initiation  and  propagation  of  Laser 
Absorption  Waves  that  have  been  observed  to  form  over  opaque  solid  targets  irradi¬ 
ated  in  air  by  an  intense  laser  beam  (Refs.  4  to  9). 

1.1  SUMMARY  OF  PREVIOUS  WORK 

The  investigation  described  earlier  was  aimed  primarily  at  achieving  an  understanding 
of  the  physics  of  Laser  Absorption  Wave  initiation  over  metal  targets  and  at  developing 
theoretical  models  for  predicting  the  initiation  phenomena  (Refs.  2  and  3).  Consider¬ 
able  effort  also  was  devoted  to  examining  the  roie  played  by  nonequilibrium  phenomena 
in  affecting  the  internal  structure  and  propagation  of  fully  developed  Laser  Absorption 
Waves  in  air. 

To  describe  the  initiation  process,  a  simple  model  of  the  laser-target  interaction  was 
postulated  in  Ref.  3.  The  model  is  based  on  uniform,  bulk  vaporization  of  the  target 
surface,  as  opposed  to  nonuniform  vaporization  such  as  might  occur  in  the  neighbor¬ 
hood  of  surface  imperfections.  As  illustrated  in  Fig.  1,  the  model  consists  of: 

•  A  layer  of  vapor  that  expands  away  from  the  surface 

•  An  outer  air  shock  layer  produced  by  the  piston  effect  of  the  expanding 
vapor 

Estimates  of  the  gasdynamics  and  of  the  laser  energy  transport  i"  the  vapor  and  air 
layers  showed  that  significant  heating  of  the  gas  by  inverse  bremsstrahlung  absorption 
takes  place  first  in  the  vapor,  rather  than  in  the  shocked  air;  and  that  the  subsequent 
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formation  of  an  absorption  wavo  in  the  air  is  connected  to  this  initial  heating  of  the 
vapor.  Nonequilibrium  phenomena  were  examined,  and  found  to  play  a  relatively 
minor  role  both  in  the  vapor  heating  process  and  in  the  physics  of  wave  propagation 
over  the  laser  intensity-wavolcngth  regime  where  subsonic  Laser-Supported  Combus¬ 
tion  (LSC)  waves  have  been  observed  experimentally  (Refs.  7,  8,  9);  namely,  for  in- 

6  2 

tensities  of  the  order  of  10  W/cm  or  less  at  C02  laser  wavelength  (10.6  n). 

In  contrast,  it  was  found  that  nonequilibrium  phenomena  are  of  major  importance  in 
the  regime  where  supersonic  Laser-Supported  Detonation  (LSD)  waves  have  been  ob¬ 
served  (Refs.  4,  5,  and  6).  Recent  computer  calculations  of  LSD  wave  propagation 
appear  to  confirm  this  conclusion  (Ref.  10). 

The  described  uniform  vaporization  model  was  used  as  the  basis  for  detailed  analysis 
of  the  laser-target  interaction  (Ref.  2).  The  analysis  is  embodied  in  a  computer  code 
called  GLIT  (Gasdynamics  over  Laser-Irradiated  Targets)  that  solves  the  one- 
dimcnsional,  time-dependent  Eulerian  gasdynamics  equations  for  both  vapor  and  air 
layers.  The  solutions  account  for  heat  conduction  in  the  solid  target,  and  for  the 
transport  both  of  laser  radiation  and  of  thermal  radiation  emitted  by  the  gas .  Code 
calculations  were  performed  for  the  conditions  of  a  typical  experiment  in  which  LSC 

waves  were  initiated  over  an  aluminum  alloy  target  irradiated  by  a  CO„  laser  at  a 

2  ^ 
peak  intensity  of  2  MW /cm  .  The  calculations  showed  the  formation  of  an  incipient 

LSC  wave  in  the  air  over  the  target  on  a  time  scale  consistent  with  the  experimental 

observations.  The  wave  initiation  mechanism  in  the  calculations  was  the  following: 

Vaporized  target  material,  heated  by  inverse  bremsstrahlung  absorption  of  laser 

energy,  emits  thermal  radiation  at  frequencies  that  are  strongly  absorbed  by  cold  air. 

An  absorption  wave  forms  when  the  air  become  sufficiently  ionized  to  absorb  laser 

energy  directly. 

These  detailed  results,  together  with  additional  approximate  calculations,  indicated 

that  the  uniform  vaporization  model  is  valid  at  CO„  laser  wavelength  in  the  regime 

2  ^ 

of  intensities  of  the  order  of  1  MW /cm  or  less  and  laser  pulse  durations  of  the  order 
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ot  1  ms  or  longer.  However,  the  model  appears  inadequate  to  explain  experimental 
data  obtained  in  short-pulse  experiments  at  intensities  between  30  and  300  MW/cnA 
Within  this  intensity  regime,  prompt  initiation  of  strong  LSD  waves  was  observed  to 
occur  in  a  shorter  time  scale  than  that  required  for  the  onset  of  uniform  surface 
vaporization  (Itcfs.  4,  5,  and  (i).  Furthermore,  this  "prompt"  initiation  was  observed 
to  occur  at  localized  "spots,"  or  sites,  that  were  randomly  distributed  over  the  irrad¬ 
iated  area  of  the  target  surface  (Ref.  4).  This  suggested  that  additional  phenomena 
were  operating  in  this  higher  power  regime. 

Examination  of  the  surface  of  aluminum  sheet  stock  of  the  type  used  in  the  experimen¬ 
tal  studies,  using  a  scanning  electron  microscope,  showed  localized  randomly  distrib¬ 
uted  pits,  scratches,  protrusions,  and  inclusions  that  could  be  the  sites  for  "prompt" 
absorption  wave  initiation  (Refs.  2,  4). 

Several  mechanisms  for  production  of  ionization  at  these  sites,  from  which  an  absorp¬ 
tion  wave  could  evolve,  appear  possible.  One  is  localized  vaporization  of  inclusions 
or  lamina  which  may  have  higher  absorptivity  than  the  surrounding  metal  surface  and 
poor  thermal  contact  with  the  bulk  of  the  metal.  This  would  lead  to  a  plasma  "microjet" 
of  the  male  rial.  A  second  mechanism  is  thermionic  emission  from  'hot  spots"  associ¬ 
ated  with  surface  imperfections.  A  preliminary  evaluation  of  this  possibility  was  made 
(Ref.  3),  and  it  was  tentatively  concluded  that  low  emission  yields  and  space-charge 
limitations  would  severely  limit  the  ionization  level  that  could  be  achieved. 

The  third  mechanism  that  was  considered  was  localized  high-field  emission  fi’om 
protrusions  on  the  surface  (Ref.  2).  It  was  found  that  protrusions  or  sharp  edges  arc 
required  to  intensify  the  normal  electric  field  of  the  incident  laser  beam  sufficiently 
to  produce  a  large  electron  emission  yield  in  the  "prompt"  initiation  laser  intensity 
regime.  Furthermore,  protrusions  of  sufficient  height  distort  the  field  enough  to 
cause  the  emitted  electrons  to  be  accelerated  to  energies  adequate  to  collisionally 
ionize  the  ambient  air  molecules.  This  negates  the  space  charge  effect  and  produces 
an  air  plasma  that  subsequently  undergoes  laser-induced  breakdown  by  cascade 
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energy  and  subsequent  heating  of  the  air  by  thermal  emission  from  the  vapor  occur 
relatively  quickly. 

Thcso  results  suggested  that,  if  the  theoretically  predicted  behavior  were  true  also 
of  the  experiments,  then  it  should  be  possible  to  correlate  all  available  experimental 
data  as  a  function  of  the  fundamental  parameters  that  govern  the  heat-conduction 
process.  Classical  heat -conduction  theory  was  used  to  deduce  the  fundamental  param¬ 
eters  as  follows:  The  GLIT  code  calculations,  in  agreement  with  the  experimental 
measurements  for  the  three  cases  cited  above,  predicted  the  LSC  wave  formation  to 
occur  during  the  initial  risetime  of  the  laser  pulse;  during  this  part  of  the  pulse,  the 
incident  laser  intensity  history  can  be  represented  by  a  ramp  function  I(t)  =  it . 

Here  t  denotes  the  time  and  i  =  const,  is  the  characteristic  slope  of  the  ramp.  For 
an  opaque  solid  target  that  vaporizes  at  temperature  Ty  and  that  has  thermal  conduc¬ 
tivity  k ,  density  p ,  heat  capacity  C  ,  and  surface  absorptivity  a  ,  the  heat- 
conduction  solution  leads  to  the  following  equation  for  the  time  of  onset  of  surface 
vaporization: 


where  is  the  initial  target  temperature  prior  to  irradiation. 

It  is  seen  from  this  equation  that  the  fundamental  laser  beam  parameter  that  governs 
heat  conduction  in  the  target  and  that  affects  the  time  required  to  produce  vapor  is  the 
characteristic  slope  i  of  the  laser  pulse.  Thus,  it  was  deduced  (Ref.  1)  that  if  the 
behavior  predicted  by  the  GLIT  code  calculations  were  true  of  the  experiments,  then 
experimentally  observed  LSC  wave  initiation  times  tj  should  correlate  well  when 
plotted  as  a  function  of  the  average  slope  of  the  laser  pulse  up  to  tho  initiation  time. 
For  simplicity,  the  following  simple  mean  value  was  employed  in  attempting  such  a 
correlation: 
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(2) 


The  correlation  is  displayed  in  Fig.  2,  which  is  reproduced  from  Ref.  1.  The  figurc 
contains  all  experimental  points  for  which  a  self-consistent  set  of  data  are  available 
(Refs.  7,  8,  11)  including  shots  for  ambient  pressures  other  than  1  atm.  For  reference, 
Eq.  (1)  is  plotted  as  the  solid  line  in  Fig.  2.  The  ciose  correlation  among  experimen¬ 
tal  data  points  in  the  latter  figure  is  evident.  Furthermore,  the  data  follow  the  trend 
predicted  by  Eq.  (1),  namely,  that  the  LSC  wave  initiation  time  varies  as  <1 >~2^3  . 

This  is  very  strong  evidence  of  the  fundamental  validity  of  the  uniform  surface 
vaporization  model  and  of  the  basic  physics  of  the  LSC  wave  initiation  process  that 
arc  used  in  the  GLIT  code.  That  is,  the  generation  of  a  vapor  layer  is  a  necessary 
precondition  to  the  formation  of  an  absorption  wave;  inverse  bremsstrahlung  absorp¬ 
tion  of  laser  energy  takes  place  in  th-  initially  weakly  ionized  vapor;  thermal  radiation 
emitted  by  the  hot  vapor  then  raises  the  air  to  temperatures  at  which  it  is  ionized 
sufficiently  to  absorb  laser  energy  directly,  and  an  absorption  wave  is  formed. 


In  addition  to  the  above  studies  of  absorption  wave  initiation  physics,  the  Midterm 
Report  (Ref.  1)  presented  an  analysis  of  a  two-dimensional  (axisymmetric)  subsonic 
laser  absorption  wave  in  air.  Preliminary  results  were  obtained  that  account  for  the 
two-dimensional  geometry  of  the  LSC  wave  propagating  along  a  finite-diameter  laser 
beam,  with  a  realistic  treatment  of  the  transport  of  laser  radiation  and  of  thermal 
radiation  emitted  by  the  hot  gas  in  the  core  of  the  wave. 

1.2  OVERVIEW  OF  THIS  REPORT 


This  report  summarizes  work  performed  during  the  concluding  phase  of  the  overall 
Laser-Target  Interaction  study  described  in  the  opening  sentences  of  the  foregoing 
introduction.  In  Section  2,  wc  develop  some  improvements  to  the  correlation  of  ex¬ 
perimental  data  on  subsonic  LSC  wave  formation  (Fig.  2).  The  heat-conduction 


.1-7 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

10(11111  *1111111  *  IMtl  COM'ANT.  INC 

i  ittntitir  it  mimn  aiiciam  coiioiaiion 


LMSC-D403747 


<i> 


rl 

J  idt  =  I(  tj  )/tj 
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Tho  correlation  is  displayed  in  Fig.  2,  which  is  reproduced  from  Ref.  1.  The  figurc 
contains  all  experimental  points  for  which  a  self-consistent  set  of  data  are  available 
(Refs.  7,  8,  11)  including  shots  for  ambient  pressures  other  than  1  atm.  For  reference, 
Eq.  (1)  is  plotted  as  the  solid  line  in  Fig.  2.  The  ciose  correlation  among  experimen¬ 
tal  data  points  in  the  latter  figure  is  evident.  Furthermore,  the  data  follow  the  trend 
predicted  by  Eq.  (1),  namely,  that  the  LSC  wave  initiation  time  varies  as  <l>  2/^  . 

This  is  very  strong  evidence  of  the  fundamental  validity  of  the  uniform  surface 
vaporization  model  and  of  the  basic  physics  of  the  LSC  wave  initiation  process  that 
are  used  in  the  GLIT  code.  That  is,  the  generation  of  a  vapor  layer  is  a  necessary 
precondition  to  the  formation  of  an  absorption  wave;  inverse  bremsstrahlung  absorp¬ 
tion  of  laser  energy  takes  place  in  thv  initially  weakly  ionized  vapor;  thermal  radiation 
emitted  by  the  hot  vapor  then  raises  the  air  to  temperatures  at  which  it  is  ionized 
sufficiently  to  absorb  laser  energy  directly,  and  an  absorption  wave  is  formed. 


in  addition  to  the  above  studies  of  absorption  wave  initiation  physics,  the  Midterm 
Report  (Ref.  1)  presented  an  analysis  of  a  two-dimensional  (axisymmetric)  subsonic 
laser  absorption  wave  in  air.  Preliminary  results  were  obtained  that  account  for  the 
two-dimensional  geometry  of  the  LSC  wave  propagating  along  a  finite-diameter  laser 
beam,  with  a  realistic  treatment  of  the  transport  of  laser  radiation  and  of  thermal 
radiation  emitted  by  the  hot  gas  in  the  core  of  the  wave. 
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solution  originally  performed  for  a  ramp-function  laser  pulse  is  generalized  so  as  to 
yield  a  rigorous  solution  for  the  laser  beam  parameter  of  an  arbitrary  pulse  shape. 

This  includes  the  singular  case  where  the  pulse  shape  function  has  intcgrablc  discon¬ 
tinuities  such  as  the  initial  discontinuity  that  characterizes  a  step-function  pulse. 

Section  3  deals  with  "prompt”  absorption  wave  initiation.  Our  earlier  model  based  on 
purely  high-field  emission  from  surface  protrusions  is  extended  to  account  for  enhance¬ 
ment  of  the  electron  emission  current  by  thermionic  emission  from  laser-heated  ma¬ 
terial.  In  addition,  following  a  recent  suggestion  by  Walters  (Ref.  4  ),  we  have  re¬ 
examined  in  sonic  detail  the  model  based  upon  purely  thermionic  emission  from  a 
relatively  smooth  surface  for  which  the  field  effect  is  negligible  at  laser  intensities 
near  the  threshold  for  "prompt"  LSD  wave  formation.  It  is  shown  that,  with  this 
model,  either  of  two  widely  different  assumptions  about  the  physical  behavior  of  the 
target  material  under  very  intense  laser  irradiation  can  load  to  predictions  that  agree 
both  qualitatively  and  quantitatively  with  the  experimental  data. 


Section  4  presents  additional  results  based  on  our  earlier  formulation  of  an  analysis  of 
a  two-dimensional  subsonic  LSC  wave  propagating  along  a  finite-diameter  laser  beam. 
A  simplified  approximate  solution  is  developed  that  demonstrates  the  significant  effect 
of  radial  flow  on  the  wave-propagation  speed. 
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Section  2 

INITIATION  OF  SUBSONIC  LASER  ABSORPTION  WAVES 

As  indicated  earlier  in  subsection  1.1,  several  detailed  GLIT  code  calculations  have 
been  performed  that  are  in  close  agreement  with  experimental  data  on  subsonic  Laser 

Supported  Combustion  (LSC)  wave  initiation  over  itanium  targets  irradiated  by  a  CO,; 

2  “ 

gas  laser  at  peak  intensities  below  a  few  MW/em  (Itcf.  1).  On  the  basis  of  the  be¬ 
havior  evinced  by  these  theoretical  eaiculations,  a  large  body  of  such  experimental 
data  have  been  correlated  successfully  as  a  function  of  the  fundamental  laser  beam 
parameter  that  governs  the  thermal  response  ot  the  laser-irradiated  target  (see 
subsection  1.1).  L’r.<  er  circumstances  where  the  laser  pulse  can  be  approximated 
by  a  ramp  function,  a  rigorous  target  heat  conduction  solution  [Eq.  (1)  ]  revealed 
this  fundamental  parameter  to  lie  tho  slope  of  the  rump  pulse.  However,  tnis  solution 
*  does  not  apply  directly  to  eases  where  the  wave  was  initiated  so  late  in  the  pulse  that 

p  the  ramp  approximation  is  unjustified.  Nevertheless,  as  shown  in  Fig.  2  of  subseelion 

1.  1,  such  cases  were  found  to  correlate  reasonably  well  as  a  function  of  the  simple 
average  slope  given  by  Eq.  (2).  In  this  section,  we  develop  a  rigorous  definition  of 
the  average  slope  that  is  valid  for  an  arbitrary  laser  pulse  shape  and  apply  this  result 
to  obtain  an  improved  correlation  of  the  available  experimental  data. 

We  begin  with  the  well-known  exact  heat  conduction  solution  for  the  surface  tempera¬ 
ture  history  of  an  irradiated  solid  (Ref.  1): 


T 


T.  =  r=-—  f  ,1-^'  -  dt' 

'  yfzkpC  J  yjl  -  V 


(3) 


where  Tj  is  the  initial  temperature  and  I  (t)  is  the  incident  laser  pulse. 
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If  the  latter  pulse  is  a  continuous  and  differentiable  function  of  time,  then  Eq.  (3)  may 
on  integrated  by  parts.  The  resulting  expression  then  may  be  solved  for  the  onset  of 
vaporization  t  defined  as  the  time  at  which  the  surface  reaches  its  vaporization 
temperature  T  .  The  solution  is  identical  to  Eq.  (2)  of  subsection  1.1,  which  applies 
to  a  ramp  pulse,  except  that  the  slope  I  =  const,  of  the  ramp  is  replaced  by  a  weighted 
average  of  the  time-derivative  i  -  dl/dt  of  the  arbitrary  pulse 


«I  » 


where 


«I  >> 


4  )1/2  i  d4 


4  =  t/tv 


In  practice,  the  laser  pulse  shape  usually  can  be  approximated  accurately  by  a 
piecewise-linear  curve-fit  of  the  form 

’  0  ,  t  <  0 

I,t  ,  0  t  <  t. 


y  i.  (t.  - 1.  , )  + 1  (t 
1  =  1  ‘ 


S-.» 


t.  -it  v  t. 

J  -  1  J 


where  we  have  defined  t.  =  0  for  j  =  0  .  The  derivative  I  is  then  piecew>se  con¬ 
stant,  and  the  integral  that  appears  in  Eq.  (4b)  can  be  evaluated  readily  to  yield  the 
following  simple  analytical  expression  for  the  average  pulse  siope 


i»  -  £  i,  [d  -  l,.,)372  -  (i  -  <j)V2]  ^  ikO  -  (0) 

j-1 
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where  the  integer  k  is  determined  uniquely  by  the  condition 


4k-l  < 


1  == 


’  k 


(7a) 


or,  equivalently, 


lk  -  l  <  lv 


(7b) 


Equation  (4b)  cannot  be  used  if  the  laser  pulse  is  either  a  ship-function  in  time,  or  is 
any  other  function  that  possesses  jump  discontinuities  at  -”hich  the  derivative  I  does 
not  exist.  An  example  that  is  of  practical  interest  for  currently  planned  experiments 
(Ref.  11)  is  the  situation  in  which  the  incident  pulse  is  tailored  so  as  to  possess  a 
sharp  initial  jump  discontinuity  similar  to  that  of  a  stop  function.  The  foregoing  anal¬ 
ysis  can  be  generalized  to  apply  to  the  situation  by  separating  the  integral  that  appears 
in  Eq.  (3)  into  two  parts 


HA 

t 

C 

t 

f  I(t' )  dt' 

lim 

f  I(tf)  dt! 

,  f  I(  t'  )  dt' 

1  ~ 

e  —  0 

J  V  t  -  V 

o  T 

J  y— t- 

where  o  <  e  <  t  .  Upon  evaluating  the  integrals  and  taking  the  e  —  0  limit,  one 
obtains  an  equation  for  t  that  is  identical  to  Eq.  (4a),  except  that  the  average  pulse 
slope  is  defined  as  follows  [compare  Eq.  (4b)]: 


« I  » 


IliiU 


-t ' 

2  l  ' 


1  - 


Ol/2idC 


(») 


where  I(  0+  )  ^  t  e  [I(  t )]  represents  the  magnitude  of  the  initial  jump  discontinuity 
in  the  laser  pulse. 
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Equations  (4b)  and  (8)  represent  rigorous  definitions  of  the  average  iaser  pulse  slope 
for  a  pulse  of  arbitrary  shape.  We  have  used  these  results  to  improve  our  earlier 
correlation  (Fig.  2)  of  experimental  data  on  LSC  wave  initiation  over  titanium  targets. 
The  results  are  displayed  in  Fig.  3,  and  include  two  additional  data  points  at  large 
«I»  that  have  recently  become  available  (Ref.  11).  As  one  might  expect,  the  two 
correlations  differ  significantly  only  near  the  lower  end  of  the  <<I  »  scale,  where 
the  observed  LSC  wave  initiation  took  place  relatively  late  in  the  pulse. 


<0.01  0.1  1  10  100 


AVERAGE  LASER  PULSE  SLOPE.  <Cl » (MW/cm2 -ms) 


Fig.  3  Correlation  of  LSC  Wave  Formation  Times  as  a  Function  of 
the  Average  Laser  Pulse  Slope  Defined  in  Eq.  (4b).  MSNW 
data  on  titanium  targets  irradiated  by  a  CO?  laser  (Rtfs.  ?, 
8,  and  11) 
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Section  3 

PROMPT  INITIATION  OF  SUPERSONIC  LASER  ABSORPTION  WAVES 

3 . 1  GENERAL  DISCUSSION 

Laser-suppoited  absorption  waves  are  initiated  at  solid  opaque  target  surfaces  in  an 
ambient  atmosphere  at  laser  beam  intensities  below  the  intrinsic  electromagnetic 
breakdown  threshold  of  the  atmosphere.  Initiation  is  accomplished  by  the  emission  of 
free  electrons  from  the  surface.  These  electrons  provide  the  "seeding"  for  a  cascade 
ionization  buildup,  whicn  eventually  becomes  a  traveling  absorption  wave.  The  mecha¬ 
nism  for  the  production  of  the  initial  "seed"  electrons  varies  over  the  laser  beam  in¬ 
tensity  range  of  concern  in  high-power  laser  applications.  At  intensities  below  a  few 
megawatts/cm2,  we  have  shown  in  the  previous  section  of  this  report  that  vaporized 

target  material  appears  to  be  the  dominant  contributor  of  the  initial  free  electrons. 

2 

However,  at  intensities  above  10-MW/em  ,  other  faster  initial  production  mecha¬ 
nisms  appear  to  become  dominant.  Two  such  mechanisms  that  we  have  investigated 
are  high-field  electron-emission  from  protrusions  on  the  target  surface,  and  thermi¬ 
onic  electron-emission  from  anomalous  patches  on  the  target  surface. 

In  this  section,  we  present  the  results  of  additional  modeling  and  analysis  of  "prompt" 

n 

initiation  of  laser -supported  absorption  waves  in  the  high -Intensity  regime  (3x  10  to 
3  x  io  W/cm  ).  Previously,  we  modeled  this  phenomenon  as  a  cascade  ionization 
process  in  the  ambient  atmosphere  triggered  by  high-field  electror-emission  from  the 
target  surface  under  laser  beam  irradiation  (Ref.  1).  In  this  model,  there  is  no  ther¬ 
mal  time  lag  Involved,  because  pure  high-field  electron-emission  depends  on  the  in¬ 
stantaneous  laser  beam  intensity  and  therefore  requires  no  surface  temperature  rise. 
However,  »f  the  target  surface  temperature  does  rise,  thermal  enhancement  of  the 
high-fteld  electron -emission  will  occur.  This  extension  to  the  pure  field-t  mission 
model  has  been  incorporated  into  our  prompt  initiation  analysis.  The  results  show  that 
if  it  is  assumed  that  the  electron-  emitting  protrusions  are  thermally  isolated  from  the 
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bulk  of  the  surface  material,  so  that  their  temperature  rise  is  controlled  by  their  ther¬ 
mal  capacity  rather  than  by  heat  conduction  to  the  bulk  of  the  target  material,  the  ther¬ 
mal  enhancement  effect  plays  a  major  role  at  the  lower  end  of  the  cited  intensity  range 
near  the  LSD  wave  formation  threshold.  However,  the  high -field  effect  alone  appears 
dominant  at  the  high  end  of  the  cited  intensity  range.  These  theoretical  results  are 
highly  dependent  on  me  geometrical  configuration  of  the  surface  protrusions.  Even 
with  somewhat  extreme  assumptions  about  the  effective  sharpness  of  the  protrusions, 
the  predictions  of  this  model  are  not  entirely  satisfactory  in  terms  of  agreement  with 
available  experimental  data. 


We  have  also  reexamined  the  purely  thermionic  electron-emission  mechanism  to  see 
[  wnether  it  offers  better  theoretical-experimental  correlation  than  does  the  high-field 

electron-emission  mechanism.  An  analysis  of  a  model  for  prompt  initiation  based  on 
thermionic  electron-emission  from  the  target  surface  an  i  concomitant  cascade  ioniza¬ 
tion  in  the  ambient  atmosphere  has  been  carried  out.  In  this  model  it  is  assumed  that 
the  thermionic  electron-emission  takes  place  from  localized  '’lamina''  (or  thin  surface 
layers)  that  are  thermally  isolated  from  the  bulk  of  the  target  material  and  that  are 
heated  by  the  incident  laser  beam.  The  emitted  electrons  act  as  the  "seed"  for  cascade 
ionization  buildup  in  the  ambient  atmospheric  region  over  the  emitting  area.  The 
results  of  the  analysis  show  that  this  thermionic  emission  model,  with  suitably  chosen 
parameters,  can  lead  to  calculated  absorption  wave  formation  times  that  agree  quite 
well  with  available  (Battelle)  experimental  data.  However,  here  too  we  must  make  a 
somewhat  extreme  assumption  about  either  the  thinness  of  these  lamina  or  the  intrinsic 
*  character  of  the  heat-transfer  process  within  the  target  material  at  very  high  heating 

rates. 


3  2  HIGH-FIELD  ELECTRON-EMISSION  FN IT IATION  MODEL 

3.2.1  Model  Characteristics 


A  detailed  description  and  discussion  of  the  model  for  "prompt"  initiation  of  laser- 
supported  absorption  waves  via  purely  high- field  electron-emission  from  protrusions 
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on  the  target  surface  has  already  been  given  (Kefs.  1,2)  and  is  summarized  in  subsec¬ 
tion  1.1  of  this  report.  However,  a  brief  review  will  put  the  latest  contribution  into 
perspective.  Figure  4  illustrates  the  model  schematically.  A  localized  protrusion, 
modeled  as  a  prolate  spheroid,  stands  out  from  the  local  "average"  target  surface. 
Coherent  Interference  between  the  incident  and  reflected  laser  beams  creates  a  standing- 
wave  electric  field  that  has  a  component  normal  to  the  surface.  This  field  is  distorted 
in  the  vicinity  of  the  protrusion,  so  that  the  normal  electric  field  intensity  at  the  tip  of 
the  protrusion  is  enhanced  (by  a  factor  F)  above  the  normal  electric  field  intensity  over 
the  "average"  target  surface.  If  the  electric  field  Intensity  at  the  tip  is  sufficiently 
large  (of  the  order  of  2  x  10^  V/cm),  electrons  are  emitted  directly  from  the  tip  (high- 
field  electron-emission).  These  emitted  electrons  are  accelerated  as  they  pass  through 
the  hlgh-field  region  that  surrounds  the  protrusion  tip.  For  protrusion  heights  less  than 
approximately  one-tenth  the  wavelength  of  the  incident  radiation,  the  kinetic  energy  that 

these  electrons  can  gain  in  traversing  this  high-field  region  is  given  by  E  h  where 

•  no  ’ 

h  Is  the  protrusion  height  and  Eno  Is  the  normal  component  of  the  standing-wave  elec¬ 
tric  field  at  the  "average"  target  surface  far  from  the  protrusion.  At  sufficiently  high 

7  0 

laser  beam  intensities  (above  1  x  10  W/cm“  for  10.6-p  radiation)  and  large  field- 
enhancement  factors  (F  «  100),  the  emitted  electrons  gain  enough  energy  to  Ionize  the 
ambient  air  by  direct  impact.  This  ionization  creates  an  air  plasma  in  the  neighborhood 
of  the  protrusion.  The  size  of  this  plasmoid  is  approximately  the  "ionization  range"  of 
the  emitted  electrons  (typically  of  the  order  of  10  p  in  air  at  S.T.  P.)  The  free  electrons 
In  this  plasmoid  then  act  as  the  priming  electrons  for  an  electron  cascade  ionization 
buildup. 

In  this  process,  the  free  electrons  absorb  energy  from  the  laser  beam  by  inverse 
bremsstrahlung,  thus  heating  to  a  temperature  high  enough  so  that  they  create  additional 
lon-electron  pairs  by  Inelastic  collisions  with  neutral  molecules.  This  is  simply  a 
laser-induced  breakdown  starting  from  substantially  preionized  air,  which  proceeds 
much  faster  than  air  breakdown  starting  from  usual  initial  conditions.  Eventually, 
the  plasmoid  electron  density  becomes  high  enough  so  that  most  of  the  intercepted  laser 
energy  is  absorbed,  and  the  plasmoid  then  grows  explosively.  The  plasmoid  is  then 
the  nucleus  of  an  absorption  wave.  We  have  assumed  in  our  analysis  that  the  "incubation 
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time”  for  absorption  wave  formation  is  the  time  required  for  the  electron  density  in 
the  plasmoid  io  increase  to  the  critical  electron  density  for  the  Incident  laser 
wavelength. 

At  the  lower  end  of  the  laser  beam  intensity  range  of  concer*'  for  prompt  Initiation  of 

7  2 

absorption  waves  (about  3  x  10  W/cm  for  COr  lasers),  the  time  required  for  the 

M 

plasmoid  electron  density  to  reach  critical  value  via  purely  high -field  electron-emission 
and  concomitant  ionization  cascade  was  calculated  to  be  sufficiently  long  to  make  it 
appear  that  some  heating  of  the  surface  and  protrusion  could  take  place .  This  would 
enhance  the  electron  emission,  and  therefore  tend  to  shorten  the  incubation  time . 

In  order  to  incorporate  this  phenomenon  into  the  prompt  initiation  model,  we  have  for¬ 
mulated  an  analytic  description  of  electron  emission  from  a  surface  protrusion  under 
the  condition  of  simultaneous  existence  of  a  very  high  normal  electric  field  and  ele¬ 
vated  temperature,  and  have  developed  a  method  for  calculating  the  temperature  rise 
of  the  protrusion  due  to  laser  energy  absorption. 

The  analysis.  Involving  considerations  of  temperature -field  emission,  protrusion  heat¬ 
ing,  and  the  numerical  results,  is  presented  In  detail  below. 

3.2.2  Analysis 

Temperature-Field  Emission,  The  emission  of  electrons  from  a  heated  surface  in  the 
presence  of  a  very  high  normal  electric  field  is  designated  as  temperature-field  (T-F) 
electron  emission  (Ref.  12).  This  is  a  particularly  difficult  regime  in  which  to  obtain 
an  explicit  analytical  expression  for  the  electron-emission  current  density  which  is 
needed  in  our  model  to  calculate  the  electron  density  buildup  in  the  plasmoid  over  an 
emitting  protrusion.  Approximations  have  been  developed  based  on  approaches  from 
either  the  purely  high-field  regime  or  the  purely  thermionic  regime  (Ref.  13) .  In  the 
former  approximation,  only  a  very  small  temperature -effect  correction  can  be  validly 
applied  to  the  high-field  emission  equation.  This  correction  term  becomes  divergent 
for  temperatures  of  concern  in  our  calculations.  In  the  latter  approximation,  the 
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well-known  Richardson-Schottky  formula  for  eleotric-fteld-enhanced  thermionic  emis¬ 
sion  is  obtained.  Unfortunately,  this  formula  does  not  apply  at  the  high  electric  field 
intensities  of  concern  in  our  calculations. 

To  cover  the  regime  between  purely  high-field  emission  and  predominantly  thermionic 
emission,  we  have  developed  an  empirical  relationship  that  asymptotically  approaches 
each  limit  and  provides  a  smooth  transition  between  them.  This  equation  gives  the 
electron  emission  current  density,  J^p  (A/cm3),  due  to  simultaneous  high  electric 
field  and  elevated  temperature,  in  the  form 


where 


JTF  “  KTF  ® 


‘TF 


k^f  s  +  kf 


TF 


-1 1/2 


k‘3  +  k’3 


Kt  =  120  5  T 


Kp  =  1.54  x  10-«W 


t  (y) 


kT  = 


KLzJd. 


8.615  x  10"5  T 
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,  6.83  x  107  <f>3/2  ,  , 

kr  =  - p — - v  ( y ) 


(15) 


.  .  .  ^  ^  1  5/3 

t  ( y )  -  i  +  q  y 


(1G) 


v(  y )  =  1  -  v 


5/3 


(17) 


y  = 


-4  1  /•> 

3.79  x  10  E 

6 


(18) 


where 

T  is  the  surface  temperature  (°K) 

E  is  the  normal  electric  field  intensity  (V/cm) 

6  is  the  work  function  of  the  surface  (electron-volts) 

When  T  is  set  equal  to  zero,  we  obtain  the  Fowler-Nordheim  formula  for  high-fitld 
emission  current  density.  When  E  is  small,  Eq.  (9)  reduces  to  the  Richardson- 
Schottky  formula  for  electric-field-enhanced  thermionic  emission  current  density. 
Equation  (9)  is  used  in  the  plasmoid  ionization  buildup  equation  in  Ref.  1  instead  of 
the  peak  current  density  formula  given  there. 


Protr  '.sion  Heating.  The  temperature  rise  of  thr  emitting  protrusion  is  estimated 
from  the  results  of  the  analyses  of  two  different  .  miting  situation  models.  In  the  first 
model,  it  is  assumed  ttiat  the  protrusion  is  in  good  thermal  contact  with  the  bulk  of 
the  target  material  so  that  both  have  the  same  temperature.  The  temperature  of  the 
bulk  material  surface  is  then  calculated  by  the  heat-conduction  model  that  wc  have 
used  for  vaporization  calculations  in  the  low-intensity  regime  (see  Section  2).  In  the 
second  model,  it  is  assumed  that  the  protrusion  is  thermally  isolated  from  the  bulk 
of  the  target  material  so  that  the  protrusion  temperature  ip  determined  by  its  heat 
capacity  and  by  the  absorbed  laser  radiation.  In  the  parameter  regime  considered  in 
this  model,  it  can  be  shown  that  the  electromagnetic  skin  depth  (depth  of  laser 


3-7 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
iocwkiio  mtsiiiii  ft  imci  coumn  ,  i  *  c 


i 


* 


LMSC-D4  03747 


radiation  penetration)  is  less  than  the  smallest  dimension  of  the  assumed  protrusion. 
The  protrusion  is  therefore  opaque  to  the  laser  radiation.  Also,  the  time  for  thermal 
conduction  through  the  protrusion  is  small  compared  to  the  time  interval  of  concern  for 
initiation;  hence,  the  concept  of  isothermal  heating  of  the  entire  protrusion  volume  is 
applicable.  Under  these  conditions,  the  temperature  T  of  a  protrusion  under  laser 
irradiation  is  given  by 


T  + 

o  p  C  V 


where 


Tq  is  the  initial  temperature 

o  is  the  effective  surface  absorptivity'  over  the  radiation  interception  surface 
area  A 

V  is  the  volume 

pC  is  the  specific  volumetric  thermal  capacity 
!  is  the  laser  radiation  intensity 


The  volume  of  a  prolate  spheroidal  protrusion  is 


,r  2,2, 

V  =  x  w  b  h 

O 


where  b  is  the  base  radius  and  h  is  the  height.  The  interception  area  presented  to 
incoming  radiation  at  an  angle  c>.  to  the  surface  normal  can  be  shown  to  be 


-  rb 


/  2  2 

b  cos  9 .  +  I  h  sin  9 .  +  b 
i  \  i 


2  „ 
coc  0 


,n 


Equations  (20)  and  (21)  can  be  combined  to  give  the  area  to  volume  ratio  in  the  form 
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cos  0.  + 


.  2 


sin 


C.  +  cos 


2 


where  x  =  h/b  is  the  aspect  ratio  of  the  protrustion. 


(22) 


Numerical  Results.  Calculations  of  incubation  time  for  prompt  initiation  of  laser- 
supported  absorption  waves  have  been  made  using  temperature -field  emission  from 
protrusions  whose  temperature  is  calculated  from  either  the  conduction  or  the  ther¬ 
mal  capacity  model . 


The  quantitative  results  are  illustrated  in  Fig.  5.  The  curves  labeled  "t  "  are  the 
times  required  to  achieve  critical  electron  density  in  a  localized  plasmoid  above  a 
surface  protrusion.  This  condition  is  implicitly  assumed  to  be  equivalent  to  initiation 
of  an  absorption  wave.  The  calculations  were  made  for  an  aluminum  surface  irradi¬ 
ated  by  a  CC2  l?«er  (  wavelength)  having  the  temporal  pulse  shape  of  the 

Battelle  laser  (Ref.  4).  The  surface  is  assumed  to  possess  a  prolate  spheroidal  pro¬ 
trusion  of  height  h(h  =  1  and  0.5  p)  and  of  aspect  ratio  1G.  For  this  geometry, 
the  electric  field  is  intensified  by  the  factor  F(  F  =  100)  in  the  neighborhood  of  the 
tip  of  the  protrusion.  The  solid  curves  (a,  b,  which  overlap)  are  for  no  temperature 
riBe  (curve  a,  pure  high-field  electron  emission)  and  for  a  temperature  rise  equal  to 
the  bulk  material  surface-temperature  rise  calculated  by  the  heat-conduction  model 
that  we  have  used  successfully  in  the  LSC  wave  regime  (curve  b,  thermally  augmented 
high -field  electron  emission) .  The  dashed  t  curves  are  fcr  the  maximum  possible 
protrusion  temperature  rise  computed  under  the  assumption  that  heat  conduction  from 
the  protrusion  into  the  underlying  material  is  negligible.  Curves  a  and  b  are  indis¬ 
tinguishable  because  the  heat-conduction  model  predicts  a  very  small  temperature 
riBe;  hence,  the  thermionic  effect  is  negligible .  The  solid  curve  labeled  "tj  67"  is 
the  time  required  to  deliver  1.G7  j/cm“  to  the  surface,  and  also  represents  a  good  fit 
to  the  measured  initiation  dwell  time  data  from  ihe  Battelle  experiments  (Ref.  4).  It 
can  be  seen  that,  although  "incubation"  dwell  times  calculated  from  the  thermally 
augmented  field-emission  model  are  not  inconsistent  with  the  experimental  measure¬ 
ments,  the  general  trend  of  the  experimental  data  is  not  well  predicted  at  the  lower 
end  of  the  intensity  range. 
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3.3  THERMIONIC  ELECTRON-EMISSION  INITIATION  MODEL 

3.3.1  Model  Characteristics 

The  prompt  initiation  Incubation  time  measured  in  the  Battelle  experiments  correlates 
very  well  with  the  time  required  to  deliver  a  fixed  energy  density  to  the  target  surface 
(Ref.  4).  This  agreement  suggests  that  the  initiation  process  might  be  related  to 
target  heating  in  a  thermal  capacity  mode . 

Walters  (Ref.  4)  examined  the  plausibility  of  this  argument  by  making  a  sequence  of 
calculations  and  empirical  correlations.  First,  he  showed  that  if  thin  "lamina"  of 
target  material  were  present  on  the  surface,  and  were  thermally  isolated  from  the 
bulk  of  the  target  material,  they  could  be  heated  to  thermioric-emission  temperatures 
(2000°  to  2500°K)  in  a  fraction  of  the  measured  initiation  delay  times  if  they  were  thin 
enough.  Second,  he  showed  that  under  the  assumption  of  steady-state  space-charge 
election  distribution,  the  electron  density  in  the  air  over  the  lamina  could  reach  values 
of  the  order  of  10  electrons/cm  for  the  thermionic  electron-emission  current 
densities  prevailing  under  the  temperatures  previously  calculated.  Finally,  it  was 
inferred  that  the  air  breakdown  time,  lor  preionized  air  with  the  above  mentioned 
"priming"  electron  density,  would  also  be  less  than  the  absorption  wave  initiation 
time  delay.  Thus,  Walters  found  no  compelling  reason  to  reject  this  initiation 
mechanism. 

We  have  reexamined  the  purely  thermionic  electron-emission  mechanism  in  order  to 
formulate  a  complete,  self-consistent  model  with  which  to  calculate  the  absorption 
wave  initiation  incubation  time  and  to  see  whether  it  correlates  with  experimental  data. 
An  analysis  of  this  model  for  prompt  initiation  based  on  thermionic  electron-emission 
from  the  target  surface  and  concomitant  cascade  ionization  in  the  ambient  atmosphere 
has  been  carried  out.  Figure  6  Illustrates  the  model  schematically.  It  is  assumed 
that  thermionic  electron-emission  takes  place  from  a  localized  thin  surface  layer  tha. 
is  thermally  isolated  from  the  bulk  of  the  target  material  and  that  is  heated  by  absorp¬ 
tion  of  a  fraction  (  a  )  of  the  incident  laser  beam  Intensity.  The  emitted  electrons  act 


LMSC-D403747 


REFLECTED 

LASER 

BEAM  / 


ATMOSPHERIC 

PLASMA 


EMITTED 

ELECTRONS 


"AVERAGE” 

TARGET 

SURFACE 


LOCALIZED 

LAMLNA 


INCIDENT 

LASER 

BEAM 


Fig.  6  Localized  Plasma  Formation  Over  Surface  Lamina 
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as  a  "seed"  for  cascade  ionization  buildup  in  the  ambient  atmospheric  region  over  the 
emitting  area.  The  absorption  wave  Initiation  time  was  chosen  to  be  that  time  at  which 
the  laser  beam  intensity  at  the  surface  decreased  significantly  from  the  incident  Inten¬ 
sity  due  to  absorption  in  the  plasma  layer  over  the  target  surface.  This  corresponds 
directly  to  one  of  the  diagnostic  measurements  used  in  the  Battelle  experiments,  where 
"significant"  absorption  was  taken  to  be  a  2-percent  decrease  in  intensity. 

3.3.2  Attenuation  Through  Plasma  Layer 

The  laser  beam  intensity  at  the  target  surface  ( I^) ,  which  is  assumed  to  be  covered 
by  a  plasma  layer,  is  calculated  by  the  WKBJ  approximation  (Ref.  2)  to  be 


\  ~  !h exp 


1 

chi  cos  6. 
ec  l 


00 

I  N  v  dz 
J  e  c 

z  =0 


(23) 


where 


*b 
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N  and  v 
e  c 
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is  the  incident  laser  beam  intensity 
is  the  speed  of  light 

is  the  critical  electron  density  for  the  incident  radiation 
wavelength 

is  the  angle  of  incidence  of  the  beam  relative  to  the  surface  i  mal 

are  the  electron  density  and  electron-neutral  collision  frequency 
in  the  plaBma  layer 

represents  the  coordinate  dimension  normal  to  the  surface 


The  fractional  decrease  (D)  in  intensity  at  the  target  due  to  dissipation  in  the  plasma 
layer  is  (1^-^)/^,  which  is 


D  =  1  -  exp 
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This  is  the  quantity  that  we  must  calculate  as  a  function  of  time  as  the  electron  density 
Increases  because  of  both  thermionic  emission  and  cascade  multiplication. 


The  power  absorbed  by  the  plasma  layer,  per  unit  area  parallel  to  the  target  surface, 
from  the  incoming  laser  beam,  is  given  by  ( -  It)  cos  0^  .  Tht3  can  be  expressed  as 
DI^  cos  0,  The  intensity  of  the  reflected  laser  beam  at  the  target  surface  is  (1  -  o  )  It, 
where  a  is  the  surface  absorptivity  of  the  emitting  area.  Beyond  the  plasma  layer, 
the  reflected  beam  intensity,  calculated  by  the  WKBJ  approximation,  is  given  by 
(1  -  0)^(1  -  D).  Therefore,  the  power  absorbed  by  the  plasma  layer,  per  unit  area 
parallel  to  the  target  surface,  from  the  outgoing  reflected  laser  beam,  is  given  by 

J(l-Of)It-(l-o)lt(l  -  D )  j  cos  =  ( 1  -  or )  ^  D  cos  flj 


Consequently,  the  total  power  absorption  from  both  the  incident  and  reflected  beams,  is 
given  by 


D  1^  cos  0j  +  ( 1  -  Qf )  Ij  D  cos 


which  can  be  written  as 


C 


DIbcos01  [l  +  (1  -<*>  <1  -  D)j 

This  power  goes  into  plasma  heating  (kinetic  energy),  vibration,  dissociation,  excitation, 
ionization,  radiation,  etc. 

3.3.3  Ionization  Buildup 

In  modeling  the  buildup  of  electron  density  in  the  plasma  layer,  we  have  made  a  number 
of  assumptions  such  that  we  overestimate  the  effectiveness  of  thermionic  initiation.  The 
electron  density  calculated  from  the  analysis  represents  an  upper  bound  on  the  actual 
electron  density  that  could  be  achieved.  This  has  the  effect  of  making  the  incubation  time 
calculated  under  these  assumptions  represent  a  lower  bound  to  the  actual  time. 
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Specifically,  we  assume  that  all  energy  absorbed  by  the  plasma  layer  from  the  incident 
and  reflected  laser  beams  is  effective  in  producing  additional  ionization.  We  neglect 
the  potential  effect  of  space-charge  in  suppressing  the  thermionic  emission  from  the 
surface,  and  do  not  consider  the  loss  of  electrons  by  diffusion  back  to  the  surface, 
recombination,  attachment  to  oxygen,  etc. 


The  number  of  electrons  in  the  plasma  layer,  per  unit  area,  due  to  thermionic  electron- 
emission  from  the  surface,  is  given  by 


t 

qe  -  I  JTdt 

t-  0 


(25) 


where  is  the  thermionic  electron  emission  rate  (per  unit  area).  This  is  just  the 
special  case  of  Eq.  (9)  for  zero  normal  electric  field.  Specifically,  this  is 

J  =  7.52  x  1020  T2  exp  j  -  - £ - r— 1  (26j 

L  8.615  x  io-0  Tj 

2 

where  J^,  above  gives  electrons/sec-cm  ;  the  other  quantities  have  already  been 
defined. 

The  number  of  electrons  created  in  the  plasma  layer,  per  unit  area,  due  to  cascade 
ionization  driven  by  absorbed  laser  power,  is  given  by  the  energy  absorbed  divided  by 
the  energy  required  to  create  an  ion-electron  pair  (  ;  ionization  potential).  This  is 


Ql  « 


_l_ 

<t>t 


f  DI^  cos  0.  +  ( 1  -  or)  (  1  -  D)  J 

t  =  0 


dt 


(27) 


The  total  number  of  electrons  in  the  plasma  layer,  per  unit  area  parallel  to  the  surface 
at  any  time  (t)  ,  is  (Q£+Qj).  We  note  that 

o 

n  .► 
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3.3.4  Surface  Heating 

To  calculate  the  thermionic  electron  emission  it  is  necessary  to  know  the  temperature 
cf  the  emitting  surface.  We  have  approached  this  calculation  from  two  widely  differing 
viewpoints.  One  is  that  there  exist  on  the  material  surface  thin  lamina,  or  "flakes,  " 
of  target  material  that  are  thermally  isolated  from  the  underlying  bulk  material.  The 
other  is  that  heat  conduction  into  the  bulk  of  the  target  material  is  rate- limited,  so  that 
the  temperature  rise  due  to  laser  power  absorption  within  the  skin  depth  layer  at  very 
high  incident  power  densities  is  affected  only  by  the  thermal  capacity  of  the  material. 


For  a  thin  flake  of  target  material  (thicker  than  the  skin  depth,  however),  heating  iso- 
thermaily  under  laser  irradiation,  we  can  use  Eq.  (19)  to  give  the  surface  temperature. 
The  area- to- volume  ratio  is  d'1 ;  hence 


T 


a  cos  0  j 

T  +  _ - 

o  pCd 


t 


t  =  0 


(29) 


where  d  is  the  thickness  of  the  flake.  The  effect  of  oblique  incidence  is  included  here 
through  the  cos  0  .  factor. 

Under  the  assumption  that  heat  conduction  is  severely  rate-limited  (that  is,  does  not 
occur  at  all),  the  temperature  of  the  surface  is  given  by 

t 

2  a  cos  0  .  f 

T  =  To  +  —6—  j  't dt  <30> 

t  =  0 
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where  6  1b  the  electromagnetic  akin  depth  in  the  material  at  the  laser  wavelength.  It 
can  be  seen  that  this  case  is  equivalent  to  the  lamina  model  with  d  =  6/2 .  That  is, 
the  temperature  rise  in  this  case  is  the  same  as  that  calculated  for  an  isothermal  flake 
of  thickness  equal  to  one-half  the  skin  depth.  The  skin  depth  (6)  is  given  by  (Ref.  14) 


6  ■■= 


(31) 


where 

X  is  the  free- space  wavelength 

c  is  the  speed  of  light 

H  and  cr  are  the  inductivity  and  conductivity  of  the  material 

It  should  be  noted  that  for  electrically  conductive  materials  such  as  metals,  a  simple 
classical  relationship  between  the  surface  absorptivity  and  the  skin  depth  exists,  and 
can  be  expressed  in  the  form 

a_  ^  4rr 

6  ~  T  (32) 

Thus,  the  ratio  a/6  in  Eq.  (30)  is  independent  of  the  electrical  properties  of  the  mate¬ 
rial  and  depends  only  on  the  radiation  wavelength. 

3.3.5  Solution  of  Equations 

The  incubation  time  for  prompt  initiation  of  laser- supported  absorption  waves  in  the 
thermionic  electron-emission  model  is  taken  to  be  the  time  at  which  a  significant  decrease 
in  laser  beam  intensity  at  the  target  surface  occurs.  In  this  context,  a  "significant1' 
decrease  is  taken  to  be  a  few  percent.  Provided  that  we  therefore  consider  values  of  D 
in  Eq.  (24)  that  are  less  than  one-tenth,  we  can  approximate  the  exponential  term  by  the 
first  two  terms  of  its  power  series  expansion,  giving 


3-17 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 


LOCIMIID 


«  i  «  1 1  i  a  i 


I  M  (  I 


COMPANY. 


I  N  c 


LMSC-D403747 


D 


00 
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C  N  COS  fl. 
ec  I 


/ 

z  =  0 


N  v  dz 
e  c 


(33) 


Furthermore,  we  assume  that  the  electron  collision  frequency  is  independent  of  position 
in  the  plasma  layer,  so  that  it  can  be  removed  from  the  integral.  Then,  using  Eq.  (28) 
gives 


D 


a; 


V 

_ c 

cN  cos0. 
ec  i 


<QE  +  QI> 


(34) 


Again,  using  the  condition  that  we  seek  a  solution  only  for  small  values  of  D,  Eq.  (27) 
can  be  approximated  by 


Q, 


(2- 


i 

a  )  cos  o.  f 

—  J  D 

t  =  o 


dt 


(?5) 


Using  Eqs.  (25)  and  (35)  in  Eq.  (34)  gives 


D 


t 

c  N  cos  0~  '  I 
ec  1  [t-O 


J  j  dt  + 


(2  -  a)  cos  O. 


D  Ib  dt 


t  =  0 


(36) 


Differentiating  this  equation  with  respect  to  t ,  and  rearranging,  gives 


.  _  i'  ( 2  -  a )  v 

dD  _  _c _ ’  ^  c  T 

dt  c  N  cp.  o  c  N  cos  0  JT 
ec  l  ec  i 


(37) 


This  differential  equation  can  be  solved  to  give 
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f  V2‘Q)  f 

catwl^i'  '  /  ^ 


- 5 _ - 

^  cos  0.  a 
ec  i 


V 

f  vQ  ( 2  -  a  ] 

tloTT!^r 


dt  I  dt 


If  we  assume  that  is  substantially  independent  of  time  during  most  of  the  ionization 
buildup,  it  can  be  removed  from  under  the  integral.  The  energy  density  delivered  by 
the  laser  beam  up  to  any  given  time  [E  (t)  ]  is  given  by 


Thus,  Eq.  (38)  can  be  written 


l 

I  ■» 


v  (2  -  a) 

D  -  exp  — 0“E(t) 

ec  vi 


r  v 

I  _ £ _ 

I  c  N  cos  9. 
.  _ «  ec  i 


v  (2-  a) 

exp  -  cN  07  E(t)  dt 
ec  u 


For  small  values  of  D  ,  the  integral  of  It  in  the  heating  calculations  of  Eqs.  (29)  and 
(30)  can  be  approximated  by  the  integral  of  1^  ;  hence 


/  ’< 


dt  £  E  ( t) 
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which  is  used  in  the  calculation  of  .  This  now  gives  the  variation  of  D  with  time 
explicitly. 


In  the  Battelle  experimental  work,  absorption  wave  initiation  time  was  measured  as  the 
time  at  which  the  laser  beam  intensity  at  the  target  surface  (observed  through  a  hole  <n 
the  target)  decreased  significantly  (approximately  2  percent)  from  the  incident  beam  in¬ 
tensity,  because  of  absorption  by  the  plasma  layer  over  the  target. 


Therefore,  in  our  analysis,  we  ask:  "At  what  time  (t)  does  D  become  significant 

a 

(say,  two  percent)?"  This  then  gives  an  implicit  equation  for  t  in  the  form 

& 


c  N  cos  6. 
ec  i 


’  v  (  2  -  a  ) 

exp  ~cN  E(ta> 

ec  M 


t  =  t 


I 


t  -0 


JT  exp 


»c(2-  cx ) 

Tn  0“ 

ec  i 


E(t) 


dt  =  0.02 


(42) 


It  can  be  seen  that  the  solution  does  indeed  depend  on  the  energy  density  delivered  by  the 

laser  beam,  but  not  in  a  very  simple  way.  Equation  (42)  is  solved  numercially  to  get 

values  of  t  for  specified  input  parameters  and  laser  pulse  shapes. 

& 

3.3.6  Numerical  Results 

Figure  7  shows  the  results  of  some  specific  solutions  of  Eq.  (42)  icr  an  aluminum  target 
irradiated  by  CO^  lasers  (10.6-ji  wavelength)  having  temporal  pulse  shapes  characteristic 
of  either  the  Battelle  or  Boeing  laser  outputs.  The  Battelle  laser  pulse  intensity  history 
was  modeled  a6  an  8-ns  delay  during  which  the  intensity  is  negligible,  followed  by  a  37-ns 
linear  increase  in  intensity  up  to  the  peak  intensity  value.  The  Boeing  laser  pulse  in¬ 
tensity  history  was  modeled  as  a  1.39-^sec  linear  increase  in  intensity  up  1 1  the  peak 
intensity  value. 
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In  these  calculations  we  have  assumed  normal  incidence  of  the  laser  beam  (  0^  -  0) ,  an 
effective  air  ionization  potential  (<^)  of  15  electron-volts,  an  effective  surface  absorp¬ 
tivity  of  five  percent  (  a  =  0.05) ,  au  initial  surface  temperature  (Tq)  of  300°K  ,  and 

an  effective  surface  work  function  (<f>)  of  3.74  electron-volts.  The  critical  electron 

19  3 

density  (N  )  for  10. 6-p  wavelength  radiation  is  1.1x10  electrons/cm  ,  and 

6C  g 

the  specific  volumetric  thermal  capacity  (  pC)  of  aluminum  is  2.35  j/cm  -  °K. 

The  electron  collision  frequency  (  v  )  in  the  ambient  air,  at  one  atmosphere  pressure, 

c  11-1 
was  calculated  for  an  electron  temperature  of  2  electron-volts  (  v  =7><10  sec  )  . 

The  "free"  parameter  in  these  calculations  is  the  layer  thickness  (d) ,  for  the  lamina 

model.  The  minimum  value  of  d  is  coi  trained  to  be  greater  than  the  electromagnetic 

skin  depth  (6)  of  the  material,  in  order  for  the  layer  to  abso/b  the  fraction  c*  of  the 

incident  laser  beam  intensity.  Thus,  the  curves  labeled  d  =  6  represent  the  physically 

realistic  minimum  values  for  t  for  a  thermally  isolated  flake.  The  other  curves  show 

the  effect  of  thickening  the  layer.  Also  shown  in  Fig.  7  is  the  time  required  to  deliver 
2 

1.67  J/cm  for  each  of  the  two  pulse  shapes.  As  noted  earlier,  this  curve  represents 
well  the  experimentally  measured  absorption  wave  formation  time  (time  for  2-percent 
intensity  reduction)  in  the  Battelle  experiments  (Ref.  4).  No  formation  time  measure¬ 
ments  were  made  in  the  Boeing  experiments.  It  can  be  seen,  then,  that  this  thermionic 
emission  model,  with  suitably  chosen  parameters,  can  lead  to  calculated  absorption  wave 
formation  times  that  agree  quite  well  with  available  (Battelle)  experimental  data. 


For  the  rate-limited  heat  conduction  model,  the  results  would  be  very  similar  to  those 
obtained  above,  except  that  there  is  no  "free"  parameter.  The  calculated  initiation  time 
would  be  the  same  as  a  d  =  6/2  curve.  That  1b,  the  predictions  for  t  for  the  rate- 
limited  heat  conduction  model  would  lie  Just  below  the  lowest  curve  (the  d  =  6  curve) 
shown  on  Fig.  7  for  each  of  the  two  pulse  shapes.  Therefore,  this  model  is  not  incon¬ 
sistent  with  the  experimental  results  when  one  considers  that  the  predicted  ta  values 
represent  lower  limits  to  the  initiation  dwell  time. 

i 

The  actual  time  required  for  prompt  initiation  of  laser- supported  absorption  waves  via 
the  thermionic  electron-emission  mechanism  will  be  greater  than  the  t  values  given 

A 

by  Eq.  (42)  for  several  reasone.  Formation  of  a  space-charge  layer  over  the  emitting 
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surface  will  suppress  the  electron  emission.  Loss  of  electrons  by  diffusion,  recombin¬ 
ation,  and  attachment  will  slow  down  the  free-electron  density  rate  of  increase.  Elec¬ 
tron  energy  loss  to  processes  other  than  ionization,  such  as  kinetic  energy,  vibrational 
and  electronic  excitation,  dissociation,  and  radiation,  will  also  depress  the  rate  of  in¬ 
crease  of  free  electrons.  Also,  the  effective  collisional  coupling  between  free  electrons 
and  vibrational  modes  in  air  may  tend  to  keep  the  electron  temperature  from  rising  as 
high  as  assumed  in  the  calculations  presented  here.  This  would  make  the  electron  col¬ 
lision  frequency  smaller.  This  would,  in  turn,  decrease  the  rate  of  energy  absorption 
from  the  laser  beam  and  consequently  slow  down  the  ionization  cascade.  In  the  rate- 
limited  thermal  conduction  model,  the  temperature  rise  of  the  surface  should  actually 
start  out  the  same  as  in  the  bulk  conduction  model  during  the  initial  low- intensity  part 
of  the  laser  pulse  risetime,  and  then  later  at  higher  intensity  become  predominantly 
heat-capacity-controlled.  This  would  lengthen  the  time  required  for  the  surface  to 
reach  the  temperature  at  which  thermionic  emission  could  begin  and  trigger  the  ioniza¬ 
tion  cascade. 

3.4  CONC LUSIONS  AND  RECOMMENDATIONS 

We  have  investigated  two  mechanisms  that  could  be  responsible  for  prompt  initiation  of 

2 

laser- supported  absorption  waves  in  air  at  beam  intensities  above  10  MW /cm  .  These 
are  high-field  electron-emission  from  protrusions  on  the  target  surface,  and  thermionic 
electron-emission  from  anomalous  patches  on  the  target  surface.  Either  of  these 
processes  potentially  could  be  responsible  for  the  production  of  triggering  electrons  that 
start  an  ionization  cascade  in  the  air  over  the  surface  under  laser  irradiation.  These 
localized  ,fbreakdown"  regions  then  develop  into  traveling  absorption  waves.  As  a  result 
of  our  analyses  of  models  of  these  two  initiation  mechanisms,  we  have  been  able  to  show 
that  neither  one  should  be  eliminated  from  consideration  until  further  experimental  evi¬ 
dence  is  available.  In  fact,  it  is  iikely  that  both  high- field  and  thermionic  mechanisms 
play  some  vole  in  the  initiation  proc  sses,  with  one  mechanism  dominant  over  the  other 
alternately  at  the  two  extremes  of  the  laser  beam  intensity  range  involved. 

The  high-field  mechanism  was  favored  initially  because  there  is  no  thermal  time-lag  in¬ 
volved  in  the  start  of  electron  emission  from  the  target  surface.  That  this  held  promise 
for  extremely  rapid  initiation  of  absorption  waves  was  indicated  by  experimental  data. 
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Although  calculations  of  prompt  absorption  wave  incubation  time  based  on  the  thermally 
augmented  high-field  electron-emission  model  are  not  inconsistent  with  the  experimen¬ 
tal  measurements,  the  general  trend  of  the  experimental  data  is  not  well  predicted  at 
the  lower  end  of  the  intensity  range.  It  should  be  pointed  out,  however,  that  near  this 
threshold  the  numerical  results  of  this  model  will  be  very  sensitive  to  the  detailed 
temporal  behavior  of  the  laser  pulse  intensity. 

The  thermionic  mechanism  was  reexamined  in  detail  after  experimental  data  became 
available  that  showed  that  initiation  time  could  be  correlated  with  time  to  deliver  a 
constant  energy  density  to  the  target,  because  this  observed  behavior  suggested  a  heat- 
capacity  effect.  Calculations  of  prompt  absorption  wave  incubation  time  based  on  the 
thermionic  electron-emission  model,  with  suitably  chosen  parameters,  agree  quite  well 
with  available  experimental  data.  Either  of  two  widely  differing  assumptions  about  the 
detailed  behavior  of  the  intensely  irradiated  target  surface  can  lead  to  the  aforementioned 
predictions  of  the  thermionically  emitted  electrons  that  induce  the  ionization  cascade. 

One  is  the  presence  of  surface  lamina,  or  flakes,  that  arc  thermally  isolated  from  the 
underlying  material,  and  that  are  of  the  order  of  the  electromagnetic  skin  depth  in  thick¬ 
ness.  This  has  been  explicitly  postulated  in  the  calculations  made  thus  far.  The  other 
assumption  is  that  heat  conduction  into  the  bulk  of  the  target  material  is  rate-limited,  so 
that  the  temperature  rise  due  to  laser  power  absorption  within  the  skin  depth  layer  at 
very  high  incident  power  densities  is  affected  only  by  the  thermal  capacity  of  the  material. 
Classical  heat  transfer  theory,  based  on  nonsaturating  (or  non- rate-limited)  thermal  con¬ 
duction,  predicts  a  transition  from  conduction-limited  to  heat-capacity -limited  tempera¬ 
ture  rise  behavior  at  much  higher  intensities  than  we  have  considered  here.  The  possi¬ 
bility  that  such  a  transition  could  occur  at  lower  intensity  due  to  nonlinear  phenomena  has 
far-reaching  fundamental  implications. 

Each  of  the  initiation  mechanisms  examined  in  our  theoretical  study  requires  the  accep¬ 
tance  of  some  important  assumptions  about  the  characteristics  of  the  target.  For  the 
pure  high-field  electron-emission  mechanism,  there  must  be  large,  very  pointed  protru¬ 
sions  on  the  surface  that  enhance  the  normal  electric  field  intensity  by  a  factor  of  order 
100.  For  thermally  enhanced  high-field  electron-emission,  these  protrusions  must  be 
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thermally  isolated  from  the  bulk  of  the  target  material.  For  pure  thermionic  electron- 
emission,  very  thin  flakes  of  target  material  that  are  thermally  isolated  from  the  under¬ 
lying  target  material  must  be  present.  These  are  all  surface  characteristics  of  the 
target.  The  rate-limited  heat-transfer  assumption,  used  in  a  thermionic  electron- 
emission  model,  is  more  fundamental.  However,  such  behavior  need  not  be  an  intrin¬ 
sic  characteristic  of  the  target  material,  but  rather  the  result  of  subsurface  lattice 
(crystal)  staicture  defects.  The  applicability  of  any  of  these  assumptions  must  be 
clarified  by  defiled  target  preparation  and/or  examination  in  all  future  experimentation. 

We  suggest  that  additional  experimental  evidence  be  gathered  to  provide  the  basis  for 
the  selection  of  the  dominant  initiation  mechanism  under  various  circumstances.  Some 
detailed  recommendations  for  such  experimental  work  have  already  been  given  (Ref.  2). 

A  critical  difference  between  the  high-field  and  thermionic  emission  mechanisms  is  the 
kinetic  energy  of  the  emitted  electrons.  If  high-field  emission  is  active,  the  emitted 
electrons  will  be  accelerated  to  energies  of  the  order  of  10  electron-volts  and  above. 

On  the  other  hand,  thermionic  emission  will  produce  emitted  electrons  with  average 
energies  well  below  one  electron-volt.  Thus,  measurement  of  the  energy  distribution 
of  emitted  electrons  (in  a  vacuum)  should  easily  differentiate  between  the  two  mechanisms. 

The  time  history  of  the  emission  current,  relative  to  the  time  history  of  laser  beam 
intensity,  will  also  help  differentiate  between  the  two  mechanisms.  As  pointed  out 
earlier,  high-field  electron-emission  essentially  follows  the  instantaneous  radiation 
intensity,  while  thermionic  electron-emission  suffers  a  time-lag  due  to  the  temperature 
response  of  the  target. 

Both  emission  mechanisms  display  a  sensitivity  to  the  angle  of  incidence  of  the  laser 
radiation.  As  the  incident  laser  beam  approaches  grazing  incidence,  the  temperature 
rise  of  the  surface  will  become  slower,  thus  delaying  thermionic  emission. 


3-25 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

LOCKMffO  ■  I  turn  *  IMd  COMPAMT.  IN( 

a  iumoiur  of  i  u  c  .  «  1 1  e  a  net  am  coiromiQN 


LMSC-D403747 


In  all  experiments  involved  in  gathering  evidence  lo  elucidate  the  nature  of  the  electron- 
emission  from  the  target  surface,  detailed  surface  conditions  and  configurations  should 
be  determined.  The  spatial  resolution  scale  must  be  less  than  1  p  (for  10. 6-pra 
radiation  experiments);  this  determination  indicates  the  extensive  use  of  scanning  elec¬ 
tron  beam  microscopy.  More  extensive  samplings  of  target  materials  and  surface  con¬ 
figurations  than  have  been  available  to  date  are  also  needed. 
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Section  4 

EFFECT  OF  RADIAL  FLOW  ON  THE  SPEED  OF 
AN  AXISYMMETRIC  SUBSONIC  LSC  WAVE 


In  the  Midterm  Term  Technical  Report  (Ref.  1),  an  analysis  of  the  axtsymmetric,  sub¬ 
sonic  Laser-Supported  Combustion  (LSC)  wave,  freely  propagating  through  air  along  a 
finite-diameter  laser  beam  was  made.  A  complete  set  of  integral  equations  which 
govern  the  two-dimensional  wave  motion  was  derived  based  on  a  generalized  integral 
method.  Results  were  presented  for  a  case  in  which  the  radial  gas  motion  was  neglected, 
but  which  included  radial  transport  of  thermal  radiation  and  heat  conduction.  This  re¬ 
sult  represents  the  solution  to  the  energy  conservation  equation  when  the  latter  is 
decoupled  from  the  radial  momentum  equation  by  the  assumptions  of  zero  radial  flow 
(pv  =  o  )  and  constant  axial  mass  flow  (pu  =  pQ  uq).  The  calculated  wave  speed  is 
therefore  that  of  a  quasi-one-dimensional  wave  and  is  about  a  factor  of  three  lower 
than  the  experimental  value. 


To  obtain  a  fully  coupled  two-dimensional  wave  solution,  the  set  of  integral  equations 
presented  in  Ref.  1  has  to  be  solved  simultaneously .  To  do  so,  we  have  followed  a 
somewhat  complicated  iteration  scheme  which  was  described  in  detail  in  Ref.  1 .  Un¬ 
fortunately,  this  multilevel  Iteration  is  very  unstable,  especially  In  the  radial  flow 
(F2  0)  iteration.  Various  stabilization  methods  were  tried,  but  without  much  success. 

In  an  attempt  to  solve  the  radial  flow  field  (F2  Q)  and  the  axial  flow  field  (F^^  Q)  without 
the  iteration  on  Fj  ^ ,  we  have  converted  the  integral  form  of  the  ^-integrated  conti¬ 
nuity  and  momentum  equation  back  into  the  differential  form.  They  can  be  written  as 


6  =  C  (l  ~  Fi,o) 


1/6 


(43) 


d  F1  0  2  F2  0 

— -r— -  0 

dx  <5 


(44) 
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"  4  7  F5, 0  +  3  (1  -  F  ~5  |16  Fl,0  F5,0  "  11  (F5,0  +  Fl,0*  "  54) 

- L0 - - - 


d  F5  0 

(42  F1  Q  +  13) - ^ 


dx 


A.  = 


1  [52QF1.0F5.0  + 

i 


86  <F5.0  *  FI,0>  +  85 


126 


F„  = 


48  Fl,0  +  12Fi.o  +  j 
63 


+  42  (1  -  F1  ^  Q)  (8  Fl,  0  F5,  0  "  F5,  0  "  2  Fl,  0  _5) 
H  F 

.  F  V/l2  5»-° 

l1  1.0/  dx 


Equations  (43)  and  (44)  are  the  first  two  ^-moments  of  the  continuity  equations  and 

Eqs.  (45)  and  (46)  are,  respectively,  the  first  moments  of  the  radial  and  axial 

momentum  equations  The  details  of  these  rj-moments  and  the  definitions  of  symbols 

were  described  in  Ref.  1.  For  a  given  F,.  Q  (X)  and  wave  speed.  Eqs.  (44)  to  (46) 

can  be  integrated  from  wavefront  (X  =  -1)  provided  that  the  initial  slopes 
dF  HV  dF  n 

_ £lS  SI  _ Li_S  and  the  constant  C  can  be  determined. 

dx  ’  dx  ’  dx 
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The  equations  are  singular  at  the  wavefront  because  of  the  following  boundary  conditions 
Fj  Q  =  1 '  F  q  =  1 ,  Fg  q  =  0 ,  at  X  =  -1  .  On  the  basis  of  physical  argument,  it 

dF5o’ 

can  be  shown  that  1 —  must  be  finite  at  X  =  -1  .  A  local  solution  to  Eq.  (45)  then 
can  be  found  near  the  wavefront.  This  solution  has  the  following  form 

1  /c 

F2,0  ~  b<'  -  Fl,0>  <«> 

Using  Fn.  (47)  ,  the  initial  slope  of  F^  Q,  Y  and  F2  Q  can  be  determined  within  a 
multiplicative  constant  b.  This  constant  and  the  constant  c  are  to  be  determined  by  the 
two  downstream  boundary  conditions,  Fy  ^  =  0  and  Y  =  0  at  X  -•  *  .  This  solution 
scheme  relieves  us  from  iteration  on  F2  0,  but  we  need  to  iterate  on  the  two  constants 
b  and  c.  In  fact,  iterations  cannot  be  avoided  because  of  the  split  boundary  conditions 
of  the  governing  equations.  In  view  of  the  extreme  instability  of  the  integral  iteration 
(on  F^  q)  mentioned  above,  this  new  scheme  seems  to  provide  a  viable  alternative. 
However,  the  limited  scope  of  the  present  study  did  not  allow  modification  of  the  com¬ 
puter  code  to  obtain  detailed  numerical  solutions  with  this  new  scheme. 

Thus,  we  have  not  been  successful  in  obtaining  a  complete  two-dimensional  LSC  wave 
solution  that  includes  the  effect  of  radial  flow.  However,  a  sufficient  set  of  trial 
calculations  were  made  with  the  original  integral  iteration  scheme  to  show  that  the  axial 
mass  flow,  normalized  by  the  free -stream  mass  flow,  is  always  less  than  unity  within 
the  wave.  Hence,  in  the  two-dimensional  case,  the  total  mass  flow  that  is  heated  upon 
entering  the  wave  is  always  less  than  that  predicted  by  the  purely  one -dimensional 
model.  This  means  that  the  radial  flow  always  tends  to  increase  the  wave  speed  over 
that  predicted  by  the  purely  one -dimensional  model.  In  the  following  paragraphs,  we 
estimate  quantitatively  the  magnitude  of  this  effect. 

Two  papers  dealing  with  the  one-dimensional  LSC  wave  appeared  recently  (Refs.  15 
and  16).  Aside  from  their  differences  in  treating  the  radiative  transport,  they  also 
differ  in  treating  the  radial  energy  loss  due  to  heat  conduction.  In  Ref.  15,  radial  con¬ 
ductive  loss  is  neglected  altogether;  in  Ref  16,  on  the  other  hand,  it  is  assumed 
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linearly  proportional  to  temperature.  We  shall  first  demonstrate  that  the  energy  equa¬ 
tion  employed  in  Ref.  15  corresponds  to  the  17-integrated  axisymmetric  energy  equa¬ 
tion,  while  the  energy  equation  employed  in  Ref.  16  corresponds  to  the  axisymmetric 
energy  equation  evaluated  at  the  axis  of  symmetry.  We  shall  then  estimate  the  effect 
of  radial  flow  in  these  two  cases  separately. 

The  axisymmetric  energy  equation  in  dimensional  form  is  (Ref.  1) 


Pu 


3h  3h  0  / k  3h\  1  3  /  k  3h\ 

3x  pv  3r  9x  ( Cp  3x  )  r  3r  y  Cp  3r  ) 


+  S  H  (r  -  Rq)  +  V  ■  Q  (48) 


After  multiplying  Eq.  (  48)  by  r,  integrating  it  from  0  to  employing  the  continuity 
equation,  and  assuming  that  r  =  0,  rpv)r  _  „  =  0 ,  —  =  constant,  we 

obtain 


^Jpu/rhdrj  = 


3x 


3  f  r  h  dr 
k  jo _ 

c  3x 
P 


OD  OO 

+  /  S  H  (r  -  Rq)  r  dr  +  /  r  V  •  Q  dr 


(49) 


where 


pu  = 


f  pu  r  h  dr 

o _ 

00 

j  r  h  dr 


In  the  one -dimensional  model,  pu  =  pQ  uq  =  constant.  In  the  axisymmetric  case, 
there  is  another  uniform  state  that  corresponds  to  the  location  where  the  final  tem¬ 
perature  is  reached  far  behind  the  wavefront  and  where  pu  =  (pu)^  =  constant.  This 
latter  constant  represents  the  mass  flow  that  enters  the  high-temperature  core  and  is 
the  dominant  parameter  that  influences  the  wave  speed.  For  estimating  the  radial 
flow  effect  on  the  wave  speed,  we  assume 
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.  r  oo  I  oo 

^  pu/  r  h  dr  I  “  (pu)^  /  r  h  dr 

*■  o  J  o 


where  (pu)  denotes  the  average  mass  flow  at  T  =  T  (X  —  *>).  Hence,  If  we  write 


WJ  UU  , 

h  =  j  r  h  dr,  T  =  J  r  —  dr 


J  c 

0  P 


K  S  =  /  r  K  S  H  (r  -  R  )  d  r  and  <fi  =  -/  rV  •  Q  dr 


Equation  (49)  is  then  identical  to  that  treated  in  Ref.  15.  Therefore, 


u  (radial  flow  included) 
o 


(one -dimensional) 


In  order  to  determine  (pu)^ ,  the  following  simplifications  are  made: 

(1)  The  wave  shape  functions  6  is  assumed  to  be  a  constant  6^  .  It  can  be 
shown  from  Eq.  (43)  that  the  variation  of  6  is  slow  except  near  the  wave 
front;  therefore,  the  assumption  5  k  ^  =  constant  far  behind  the  front  is 
reasonable . 

(2)  The  functional  relationship  between  the  pressure  field  and  the  axial  flow 
momentum  is  given  by  the  one-dimensional  axial  momentum  equation. 

(3)  pv  is  assumed  to  be  small  so  that  it  is  neglected  in  the  radial  momentum 
equation . 

Under  these  approximations,  the  continuity  and  momentum  equations  in  dimensional 
form  are 


+  (vpv)  -  o 

dx  <5  ti  drj  '  ’ 

00 
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where  rj  =  y 
00 


pu2  +  P  =  pQ  u2  +  pQ 


9v  1  0p 

pu^- =  r  b? 

oo  ' 


(53) 


(54) 


By  combining  Eqs.  (  53)  and  (54)  ,  we  obtain 


&v  J_  f  I  fy>u  ,  9ul 

8x  50olP  ^ 


(55) 


Equations  (52)  and  (55)  are  to  be  solved  for  pu  and  pv  as  functions  of  p  ,  subject  to 
tlie  following  boundary  conditions;  at  x  —  -  «  ,  pu  ^  p  u  ,  pv  =  0  ;  at  x  -*  «  ,  v  -  0, 

dx  dx  U* 

As  outlined  in  Ref.  1,  we  assume  the  radial  profiles 


pu  =  p. 


=  0„  %  "  <’  -  l/»  1  >  F2,  0 


"o(Pl.O  +  <l  '  Fi,o>"2) 

(\o  +  <'  -F4.0 ?=^(F5,0  +  '1  -F5.0>’>2) 


Based  on  the  profiles,  we  approximate  pu  as 


(56) 


}  V  P  u 


hd  r; 


pu 


_  iMF1.0F4.0+l/!!<ri.0+F4.<>>  +  1l<’oUo 


/ 


1  +  F 


rj  h  d  tj 


4,0 
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Since  F4  represents  the  dimer sionless  enthalpy  at  the  line  of  symmetry,  it  is 
directly  proportional  to  the  temperature.  From  the  one-dimensional  solution,  we  have 
Tb  »  Tq.  Therefore.  ^F^  ^  »  1,  and  we  obtain 


(P^.0 


2/3 


[(F>.»L + 1/2  ] 


P  U 
Ko  o 


(57) 


where  denotes  Fj  Q  at  x  —  *> 

We  note  that  at  x  ■—  °°  ,  v  —  0,  and  0  for  all  17;  hence,  to  find  ,  it  is  only 

necessary  to  set  the  rj-moment  of  the  right-hand  side  of  Eq.  (55)  to  zero,  i.e., 


Thus,  we  obtain 


Upon  assuming  a  perfect  gas  equation  of  state,  we  have  F5  q  ~  F4  0;  hence- 
(F5,o)„  »  >•  sothat 

(F  o)  “  l/«  (58) 

p  U  q 

Upon  substituting  (58)  into  (57)  ,  we  find  The  effect  of  radial  flow  in 

this  case  is  therefore  to  increase  the  wave  speed  as  given  in  Ref.  15  by  a  factor  of  y/4. 
Comparisons  between  wave  speeds  are  shown  in  Fig.  8. 
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Fig.  8  Propagation  Velocity  Versus  Laser 
Flux  Corresponding  to  a  C09  Laser 

It  is  of  interest  to  obtain  an  approximate  solution  for  the  axial  mass  flow  Q  and 
the  r-derivative  of  the  radial  flow  F9  „  on  the  line  of  symmetry.  This  solution  will 

fc.  ,  v 

demonstrate  the  significant  difference  between  the  one-  and  two-dimensional  cases. 


After  substituting  Eq.  (56)  into  Eqs.  (52)  and  (55),  multiplying  the  resulting  equation 
by  T]  ,  and  then  integrating  them  from  p  =  0  to  p  =  1,  we  obtain  the  governing 
equations  for  Fj  Q  and  Q  as 


d  F 


dx 


1J) 


3  <5 


2,0 


(59) 
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Fig.  9  One -dimensional  Subsonic  LSC  Wave  Temperature  Profile 

)  To  obtain  an  analytic  solution  to  Eqs.  (59)  and  (60)  ,  the  temperature  profile  is 

assumed  to  be  approximated  by  the  dotted  line  as  shown  in  Fig.  9.  Now,  for  con¬ 
stant  pressure,  Fj.  _  ~  T,  we  therefore  assume 


F5,0  1 


x  <  0 


F5,0  (F5,o)  P>  1  x  >  0 


Under  this  assumption  and  the  conditions  that,  at  x  —  -  Fj  Q  =  1  ,  — ^  -  =  0  ,  at 

dFl  0  *  dF10<0+>  dF,  n(°') 

x-.  =  0,  atx  =  0  F,  0(0*)  =  F10(0-)  and - =  ^ - 

obtain  the  following  solutions  for  Fj  Q  and  F^  Q 
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F1,0  '  1  -  °-452  e 


J_  /so . 
dJ  3 


=  0.382  e 


+  1/6 


x  <  0 


x  >  0 


0. 619  e 


=  0.619 


x  <  0 


x  >  0 


These  solutions  are  plotted  in  Fig.  10  for  6^  =  l;  we  see  that  the  axial  mass  flow 
is  substantially  reduced  due  to  radial  flow. 


Fig.  10  Axial  Mass  Flow  (Fjq)  and  r-Derivative  of  Radial  Flow  (F2o) 
on  the  Line  of  Symmetry 
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We  now  examine  the  behavior  at  the  line  of  symmetry,  wherv  pv  -  0  ,  Equation  (  4b) 
reduces  to 


ah  a  /.  3T\  _  k  a^h  ...  dg  ^dQ 

=  •s(ks)t  2^Z2(0'  +  as  +  s  ^ 

p  9r 
k  a2h 

If  we  write  p  u  ~  (pu)^  =  constant,  2  - - 1  (0)  -  -j  ,  S  =  <(>  and  Q  =  -F  , 

p  an 

Eq.  (62)  becomes  identical  to  that  treated  in  Ref.  16.  The  wave  speed  which  includes 
the  radial  flow  effect  is  again  given  by  Eq.  (  51 )  . 

In  this  case  (pu)^  is  to  be  determined  by  setting  the  right-hand  side  of  Eq.  (55)  to 
zero  at  line  of  symmetry,  i.e., 


1  fl  apu  9ul 

0?8"  aA=o, 


=  0 


(63) 


X— 


Because  of  the  ellipticity  of  the  governing  equations,  profiles  in  t)  are  required  even 
though  only  the  line  of  symmetry  solution  is  sought.  We  again  assume  the  profile 
given  by  Eq.  (56).  From  (56)  ar.d  <6?  we  obtain 


*  "o  %  (Fl,  o)  '  %  uo 


(M. 


(F5,c)  ' 


Again,  since  Fr  Q  »  1 ,  we  have 

J  ’  CO 


P0U0 


(64) 


% 


»M(  | 
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The  one-dimensional  wave  speed  given  in  Ref.  15  is  to  be  increased  by  a  factor  of  3/2 
in  order  to  account  for  the  radial  flow  effect.  Comparisons  between  wavespeeds  are 
listed  in  Table  1  for  those  cases  treated  numerically  in  Ref.  13. 


Table  1 

COMPARISONS  OF  THE  TWO-DIMENSIONAL  WAVESPEED  WITH  THE 
RESULTS  FROM  THE  ONE-DIMENSIONAL  MODEL  OF  REF.  16 


Case  1 

Uo 

Case  2 

Case  3 

One-dimensional,  linear  theory  (Ref.  16) 

9.36 

62.5 

35 . 75 

Radial  effect  included  (Eq.  (64)  1 

14 

94.4 

50.6 

Ref.  16 

40 

75 

57 

In  conclusion,  from  the  comparisons  shown  in  Fig.  8  and  Table  1,  we  see  that  the 
omission  of  radial  flow  in  the  one-dimensional  LSC  wave  model  may  be  largely  respon¬ 
sible  for  the  sizable  discrepancy  between  the  one -dimensional  results  and  experi¬ 
mentally  measured  wavespeeds. 
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